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Abstract 

We study decay rates and spectra of B — > Apir, T,°pir, T,~ nn, H°E + 7r, E~ S° it and H~A7t 
modes under a factorization approach. The baryon pairs are produced through vector, axial vector, 
scalar and pseudoscalar operators. Previous predictions, including ours, are an order of magnitude 
too small compared to experiment. By incorporating QCD counting rules and studying the asymp¬ 
totic behavior, we find an earlier relation between the pseudoscalar and axial vector form factors 
to be too restrictive. Instead, the pseudoscalar and scalar form factors are related asymptotically. 
Following this approach, the measured Ap7r + rate (~ 4.0 x 10~ 6 ) and spectrum can be understood, 
and A should be dominantly left-hand polarized, while we expect B(T,°pir + ) ~ 1.6 x 10~ 6 . These 
results and other predictions can be checked soon. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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I. INTRODUCTION 


Several three-body baryonic B decays such as B —> D* p n 


ll\,PpK B and D* pp 0 have 


emerged recently, even though there is only one single two-body baryonic mode B° A+ p 


that is observed so far 
enhanced over two-body 



It has been argued that three-body baryonic modes could be 
, by reducing energy release to the baryons via emitting a fast 


recoil meson. One consequence is enhancement near baryon pair threshold in three-body 


modes. In our study of B° —» D*~ pn 0, assuming factorization , we obtained ~ 60% of 
experimental rate from the vector current contribution, and the decay spectrum exhibits 
such threshold enhancement. The same threshold enhancement effect was predicted for the 
charmless ppii mode Q, and, interestingly, the newly observed first ever charmless baryonic 


mode, B —► ppK, showed similar feature [2], The measured decay rate can be understood 


to some extent [9j and the spectrum can be reproduced by using the factorization approach 


and QCD counting rule arguments 


Q 


Other charmless modes such as ppn, Apn, T, 0 pn 


have been studied under the factorization assumption and B(Apn + ) = (3-5) x 10~" and 
B(T,°pn + ) = (0.8-1.8) x 10 -6 were predicted Jl, 10]. 

Recently, Belle reported jll| 

B(Apn + ) = (3.97^ ± 0.56) x 10~ 6 , (1) 

and B{TA pn + ) < 3.8 x 10“ 6 at the 90% confidence level. While the Apn decay spec¬ 
trum exhibits threshold enhancement as expected, the measured rate turns out to be an 


order of magnitude higher than predicted [3,|1Q[. Furthermore, previous predictions placed 
B(YPpn + ) considerably above B(Apn + ). If the factorization approach is not to be aban¬ 
doned, where could things go wrong? 

We had noted that the B —> Apn + mode is sensitive to how one treats the vacuum 


to A p pseudoscalar matrix element [l(] under factorization. The analogous situation for 
meson case is known to be enhanced. In this work, we revisit these two modes, as well 
as some SU(3) related modes such as n n, S° £ + n, 5“ S° n and S - An. With the help 
of QCD counting rules and taking into account the asymptotic behavior of baryonic form 
factors, we can now account for the observed Apn rate and spectra, where A p production is 
dominated by the pseudoscalar density. After improving the situation on the Apn rate, we 
study the A polarization, which is known to be useful for constructing CP- and T- violation 


0 


observables |6] - We are able to make some predictions as well. 


2 






Our formulation is given in the next section, followed by results and discussion. 


II. FORMALISM 


Under the factorization assumption, the three-body baryonic B decay amplitude consists 
of two parts. For one, the baryon pair is current-produced in association with a B to meson 
transition. For the other, the B makes a transition to a baryon pair and the recoil meson 
is current-produced jicj]. The B —> ppK mode receives both contributions, but the Ap^ r + 
mode, and analogously its SU(3) related modes such as E°p7T + , E - fn: + , 5° E + 7 r + , 5“ E° 7r + 
and A7 t + , receive only the current-produced contribution. We shall apply the term 
“current-produced” to scalar and pseudoscalar densities as well. 

Take, for example, the B° —> Ap7T + decay. Under factorization, the amplitude is |l^| 

M(Apn + ) = ^|(7r + |h7 /i (l - 75) & 15 °) / (V ub V* s ai - V tb V* s a A ){Ap 1^(1 - 7 5 )m|0) 


-\-2a§VtbV t 


(PA + Pp)n 


ts 


<Ap|s(l +75HO) 


( 2 ) 


m b — m 

The baryon pair Ap is produced from vacuum through 37 ^( 75 )^ and s {p/ b )u operators, 
while B° to 7r + transition is induced by U'y^b current. Note that, from isospin symmetry, 
we have (n 0 \u'y fJ- b\B~) = (7r + \u r y fl b\B 0 )/\/2, hence Ai^Apn 0 ) = M.(Apn + )/\/2. For these 
current-produced modes, we have 

r B 


B{B~ 


BB 7T°) = 


2 t b 0 


B{B (] 


BBtD 


(3) 


where r-g 0 B - are the B° and B meson lifetimes, and B stands for some baryon. 
The (axial) vector current-produced matrix elements are decomposed as 

W) 


<BB|U m | 0 > = u(j>b) UF 1 + F 2 )^ li + 

(bb'|a m |o) = u(p B ) {gA{t)'r^ + 


m B + m B 

h A ( t ) 


Pb) u f V M > 


(PB+Pg') fl 75%), 


(4) 

(5) 


m B + m B / 

where F L2 , g A and h A are the induced vector (Dirac and Pauli), axial and the induced 
pseudoscalar form factors, respectively, and t = (p B + p W ') 2 = m 2 _,. The scalar and 

4A BB 

pseudoscalar matrix elements associated with the a$ term of Eq. © are expressed as 


(BB \S\0) = f s (t) «(pbMPb')> 
(BB , |P|0) = g P (t) u(pb)75^(Pb')- 
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( 6 ) 

(7) 






It is the gp(t) form factor that is the focus of our attention, where we offer a refined discussion 
in face of B —> A p n + data. 

The scalar and vector matrix elements can be_ related by the equation of motion, 


(BB 4'V/|0) = i (rri q — m q >)(BB \qq'\0), giving 


aQ 




m n 


m q > 


( 8 ) 


rri q i as well. For 


We note that it is safe in the chiral limit m q , m q i —> 0, and for m q 
example, for (Ap|sn|0) we have (m\ — m p )/(m s — m u ) ~ 1. For the modes studied here, the 
factor (m B — rriB')/{m s — m u ) varies by 30, 40%, which illustrates SU(3) breaking. 

The pseudoscalar and axial current matrix elements can be analogously related. Using 
(BB 7 1 d^A^O) = {m q + m g /) (BB^gi y 5 q'\0), we have 


9 a (t) + 


, . . m 0 + m a / , . 

h A {t) = - s - l — q -gp(t). 


( 9 ) 


(m B +m B 0 2 7 rriB+m-B 
As m q , m q i —> 0, we get —> — g A (t ) (m B + m B /) 2 /t Q. Since the m g /m B ratio is small, 

one is close to the chiral limit, hence the dependence of hAt) on gp(t) is weak. However, to 
ensure good chiral behavior, we previously followed Ref. 9( and took 

m GB ( m B + rn w ) 


Q 


9p(t) = ~9A(t) 


( 10 ) 


(m q + m q '){t - m 2 GB Y 

where tugb is the corresponding Goldstone boson (e.g. kaon) mass. That is, gp{t) is obtained 
by changing the 1/t term in the asymptotic form of h A {t) to 1 /(t — m,n p) and make use of 
Eq. m Indeed, Eq. m gave too small a rate for B —> A p n + |9, [lfl|. Due to the small 
quark-baryon mass ratio in Eq. ©, we note that g A and h A are insensitive to gp. Therefore 
in the previous approach we need a very precise information on both g^ and h A , which is 
unavailable so far, to pinpoint gp. In this work we choose a different strategy by studying 
g A ,p directly and obtaining h A throu gh E q. @. 


According to QCD counting rules |l2j], both the vector form factor F\ and the axial form 
factor g A , supplemented by leading logs, behave as 1/t 2 in the t —► oo limit. This is because 
we need two hard gluons to impart large momentum transfer. Similarly, considering the 
bilinear structure of the S and P operators, the scalar form factor f s and pseudoscalar 
form factor gp also behave as 1/t 2 in the asymptotic limit. However, due to the need for 
helieity flip, one needs an extra 1/t for F 2 and h A , hence they behave as 1/t 3 . We see that 
Eq. (gOP gives 1/t 3 rather than 1/t 2 asymptotic behavior for gp, which is symptomatic. In the 
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TABLE I: Relations of baryon form factors F\ + F 2 , gA and gp with the nucleon magnetic form 
factors G p ^, Da,p and Fa,p via the ( su)v,a,p operators. Replacing G^ 1 by G p p n in the second 
column, one obtains F\ + F 2 1 /(wb + m^/) 2 . 


bb' 

Fi + F‘2 

9A,P 

A p 

[3 r<p 
y 2 

“Te ( D + 3F )a,p 

T,°p 

^(^+20.,) 

75 

T,~ fi 

~( g m + 2G m) 

(D - F) a ,p 

[i] 

o 

Ml 

+ 

piP pm 

U M U M 

{D + F) a ,p 

[i] 

M| 

o 

1 (pP pm \ 

75 \yM <j m 1 

75 ( d + F)a,p 

E~ A 

\/i { g m + g m ) 

-75 < D - 3F Up 


electromagnetic current case, the asymptotic form has been confirmed by many experimental 
measurements of the nucleon magnetic form factor G= Ff'" + F p,n , over a wide range 
of momentum transfers in the space-like region. The asymptotic behavior for G p , r also 
seems to hold in the time-like region, as reported by the Fermilab E760 experiment jl3| for 
8.9 GeV 2 < t < 13 GeV 2 . Another Fermilab experiment, E835, has recently reported fldl j 
G p m for momentum transfers up to ~ 14.4 GeV 2 . An empirical fit of \G P M \ = Ct ~ 2 [ln(f/Q 2 )]- 2 
is obtained, which is in agreement with the QCD counting rule. 

The current induced form factors F \, F 2 for the modes studied here can be related to 
the nucleon (Sachs) magnetic and electric form factors Gm,e, as shown in Table Q], where we 
also give the SU(3) decomposition of gA and gp in terms of the form factors Da,p and Fa,p- 
The F] + F 2 terms are in fact obtained by using 


D v = --G 


Mi 


Fy = G p M + -G 


T M1 


( 11 ) 


with SU(3) decompositions similar to that of (?a,p- We can decompose fs similarly into Ds 
and Fs, with (compare Eq. (JHJ)) 


D s 


m-B — m-B> 
m q — m q t 



me-mg. / 1 \ 

rn q - rn q i \ 2 ) 


( 12 ) 
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From the factorization assumption and Tabled we expect 

B(Ap7i + ) ~ 2 B(Ap7i°), 

£>(E“ n7r + ) ~ 2/5(E“ h7r°) ~2/5(E°p7r + ) ~ 4/5(E°p7r°), 

£>(5° E + 7 r + ) ~ 2 / 5 (H°E + 7 r 0 ) ~ 2/5(5“ E° vr + ) ~ 4/5(5“ E° vr 0 ). (13) 


There is considerable data on the nucleon magnetic form factors. This allows us to make 


a fit 


Q: 


= E w 


i =1 


1111 -k 




2—1 


1111 h 


-7 


(14) 


where 7 = 2.148, 27 = 420.96 GeV 4 , x 2 = -10485.50 GeV 6 , x 3 = 106390.97 GeV 8 , x 4 = 
-433916.61 GeV 10 , x 5 = 613780.15 GeV 12 , y { = 292.62 GeV 4 , y 2 = -735.73 GeV 6 , and 
A 0 = 0.3 GeV. They satisfy QCD counting rules and describe time-like electromagnetic 
data such as e + e“ —> NN suitably well. The data is extracted by assuming \G P E \ = \G P M \ 
and | G e \ = 0 (which gives better fit compared to the \G E \ — |G’V case Q). With the fit 
of Eq. (fUD , time-like G(y is real and positive (negative) lol Il7l|. It is interesting to note 
that the fit coefficients Xj.s alternate in sign, and likewise for r/jS. Just two terms suffice for 


to perturbative QCD [3], asymptotically (t 00 ) one expects G'h/G^ = —2/3. We find 
that the fitted parameters for G'" f with the \G%\ = 0 assumption gives G^/G/( f —> —yi/xi = 
—0.70, which is within 5% of the QCD expectation. Note that, by use of Gm = F 4 + F 2 and 
asymptotically F 2 jF\ —>■ 1/t —> 0, we have F{ l /F p —> G n M jG p M —> —2/3 as well. 

The F 2 term can be related to (Ge — Gm)/[ t/(mB + m B') 2 ~ However, we do not have 
much data on time-like nucleon Ge■ Thus, we concentrate on the F\ + F 2 term in Eq. (jij) as 


the latter because the neutron magnetic form factor data is relatively sparse [ 7 ]. According 


we did in Refs. |7l llfl|. We also use Gm in place of iq in Eqs. (JHJ) and The effect of the 
F 2 (or equivalently Ge — Gm) contribution can be estimated by using form factor models 
such as Vector Meson Dominance (VMD), where both Ge and Gm are available. 

The time-like form factors related to Da, Fa are not yet measured, but, as pointed out 


in Ref. 


• 2 ], 


their asymptotic behavior at t —> 00 are known Hi and useful. Asymptotically, 


they can be described by two form factors, depending on the reacting quark having parallel 


or anti-parallel spin with respect to the baryon spin {TdI ] . 
factors in terms of G'((;” as t —> 00 , one has 


By expressing these two form 


D a ^G p +-G 


t Mj 


Fa 


Lx » 


M 


G 


M- 


(15) 
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In similar fashion, in the asymptotic region the fs and g P form factors for the chirality flip 
operators S and P can be expressed by just one form factor, with spin of the interacting 
quark parallel to the baryon spin. Anti-parallel spin corresponds to an octet-decuplet instead 
of an octet-octet baryon pair. Since gp (equivalently D P , F P ) and fs are related to the 


n 


same form factor, by following the approach of Ref. 19|, as shown in Appendix A, we have 

D P (s) 


9p —► fs, 


3 

2 ’ 


(16) 


Fp(s) 

as t —> oo. This is a non-trivial requirement and it is not obeyed by Eq. cm We note 
that Eq. m is obtained without the use of the equation of motion. The requirement of 
D s /F s 3/2 is consistent with Eq. (H21h which follows from Eq. ((SJ) by using Fj‘/Ff —> 
G n M /G p M —■► —2/3 asymptotically 18j. Thus, the use of the equation of motion for fs in 
Eq. (0) is consistent with the asymptotic relations in Eq. cm 

The asymptotic relations hold for large t, hence they imply relations on the leading terms 
of the corresponding form factors. In general, more terms would be needed. In analogy to 
the neutron magnetic form factor case, we express Da,p, Fa,p up to the second term jlDi ]. 



1 “7 


-7 


-7 


-7 


(17) 

The asymptotic relations of Eqs. cm m imply di — xi — 3 yi/2, fi = 2aq/3 + yi/2, 
di = (3yi/2)[(m B - m B ')/(rn q ~ m q ')\ and /i = {x 1 - yi/2)[(m B - m B >)/(m q - m^)], while 
further information is needed to determine g? 2 , fi, di and c^, as we will discuss in the next 
section. We note that the anomalous dimensions of g P and fs may not be the same as that 
of F\ 2 and gA■ However, their effect is logarithmic hence not very important, and we apply 
the anomalous dimension of F\ to others for simplicity. 


It is useful to compare with Refs. QQ on the treatment of g p (t) (or equivalently on 
liA(t)), namely Eq. (HU1) . As a working assumption, this form of g P (t) with ut.q B /( m q + m q i ) 
factor was useful in particular for the good behavior of the pseudoscalar matrix element in 


7 














the chiral limit. However, it may be too restrictive in three aspects: gp(t) oc gA(t), which 
is too strong an assumption; the appearance of the Goldstone boson pole in time-like form 
factors, although itiq b is way below baryon pair threshold; and a 1 /f 3 asymptotic behavior, 
rather than the 1/t 2 form as expected from QCD counting rules. Ultimately, it does not 
satisfy the asymptotic relation of Eq. (USD- We have improved on these points in our present 
treatment of gp(t). 


III. RESULTS 


It is straightforward to use Eq. © to calculate B —> Apn and similar rates. Before we 
start, let us first specify the parameters used. We take (ps (or 7 ) = 60° Q and central 


values of \V c b\ and \V u b\ from Ref. 
and rm = 4.88 GeV at g = 2.5 GeV 
f. [23]- 


We use rn u (d)/m s = 0.029(0.053), m s = 120 MeV 
22]. The B —> 7 T transition form factor is given in 



We 

21 , 

22 


Ref. 


For effective Wilson coefficients, we use cq = 1.05, 04 x 10 4 = —387.3 — 1217 and 
a,Q x 10 4 = —555.3 — 121? from Ref. |2j| with N c = 3. 

Following Ref. [Kj, we use the axial vector contribution to 5° ^ T* pn decay to 
constrain / 2 and d 2 . Since there is no scalar and pseudoscalar contribution in this tree 
dominated mode, we simply use the chiral limit form of hj\(t) = —^(f) (■ m p + m n ) 2 /t. The 
gp contribution is suppressed by the quark-baryon mass ratio. We update our previous 
calculation 3 ] using the present input parameters, finding the vector part of the branching 
ratio to be B V {D*~ pn) = 11.9 (a® ff /0.85) 2 x 1CT 4 , where the same a® ff value as in Ref. 25] is 
used. To reach the central value of the measured rate B(B° —> D*~ pn) = (14. 5+30 ± 2.7) x 
1CT 4 [l|, using d- 2 + / 2 = —956 GeV 6 0, we find Ba(D*~ pn) = 2.6 (a® ff /0.85) 2 x 1CT 4 from 
the axial current. Although the value of <i 2 + / 2 is about half of what was used in Refs. 0] 
and 0 ], the change only affects the branching ratios of the charmless modes studied here 
at the 1CU 8 level. Following Ref. 0], we use d 2 = / 2 . 

With the axial contribution fixed, and with the scalar and vector contribution related by 
the equation of motion (Eq. we give in Table HT1 the vector plus scalar contribution (By) 

and the axial plus pseudoscalar contribution (Ba) to B° —> BB ir + branching ratios. For Ba, 
we show two cases with either vanishing or non-vanishing d 2 and / 2 from the pseudoscalar 
form factor, which is yet to be fixed. Since the contribution from the vector plus scalar part 
does not interfere with the axial plus pseudoscalar part, the branching fraction is a simple 
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TABLE II: Branching fractions for BB7r + modes arising from the vector and scalar parts (By), 
and from the axial and pseudoscalar parts ( Ba )• The latter are given for the two cases of using 
the asymptotic gp (d 2 = ]i — 0) or the fitted gp (d 2 = f 2 = — 952 GeV 6 ) from the Apn rate. The 
branching fraction is a simple sum of the two, i.e. B = By + Ba- Rates for BBV 0 modes are about 
one half of those shown. 


Modes 

£>y(10 _6 ) 

Ba{ 10- 6 ) 

use asymptotic gp 

use fitted gp 

Apir + 

0.13 

7.97 

3.84 

pTT + 

0.88 

0.70 

0.70 

E - fnT + 

1.79 

1.41 

1.41 

[ 1 ] 

0 

Ml 

+| 

+ 

0.17 

2.23 

1.20 

3“ S°7T + 

0.09 

1.14 

0.63 

[I] 

1 

>1 

+ 

0.15 

0.38 

0.20 


sum of the two, i.e. B = Bv + Ba, just as for B° —> D* pn. By using the relation of Eq. ©, 
^(BBV 0 ) can be read off from Table ITTl by a simple 1/2 factor. 

We find B v ( A[E°]h7r + ) = 0.13 [0.88] x 10 -6 . We note that B v (Apn + ) is consistent with 
previous studies |3, lOj, while By(YP pir + ) becomes slightly larger because of the different 
input values of the neutron magnetic form factor parameters (r/;s). Clearly, By(Ap7i + ) part 
is still an order of magnitude below the measured [llj branching ratio of Eq. Q. Before 
invoking the pseudoscalar form factor of Eq. (El), let us make sure that other modifications 
are insufficient for the order of magnitude difference. 

Recall that in the vector and scalar sector, we concentrated on iq + F 2 contributions 
without including the Gp — Gm effect since Gp data is unavailable. As noted earlier, one 
can try to estimate the Gp — Gm effect by using some form factor model where both Gp 
and Gm are given. We use a VMD model, Ref. [la], which was discussed in our previous 


work 


. Since iq Ap (f) and F 2 p {t) can be expressed in terms of G p hr and G P E , and since the 


? Ap 


u 


VMD model describes G P M data better than G'lj (time-like) data [16|, perhaps the Ap^ 
mode may be a better place to estimate the Gp — G]\f effect. By incorporating VMD with the 
previous section (following similar approach of Ref. |2j), we obtain By(Apn + ) = 0.27 x 10~ 6 . 
Although we gain by a factor of two compared to Table El the effect is still of order 10 7 , 
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FIG. 1: (a) dB{Ap'K + )/dm\p spectrum, where solid (dashed) line is for using the fitted (asymp¬ 
totic) gp of d 2 = I 2 = — 952 GeV 6 (0); (b) dB{TiN tt + )/ dm^fj spectra, where solid (dotted) line is 
for ( Ti~n ). The plots for 7r + replaced by 7 r° are expected to be similar but factor of 2 lower. 


and is insufficient to account for the measured Apn + rate. The effect of Ge — Gm is not 
likely to fill the gap between B\/(Ap7T + ) and the measured B(Ap7T + ). 

We thus need to turn to the axial and pseudoscalar contributions. Let us start by using 
only the d\ and fi terms of gp determined by the asymptotic relation of Eq. (ESP, i.e. taking 
d -2 — f '2 — 0. It is remarkable that, as given in the first case for Ba in Table |E] (column 
three), the 1/t 2 terms of Dp and Fp alone give B(Apn + ) ~ 8 x 10 -6 , or overshooting the 
experimental value by a factor of two! This is striking compared with the previous calculation 
using the ansatz of Eq. m, which gave results an order of magnitude too small P, [itj . 
Now, we know that the sign of x,s and 7/jS alternate hence Gm gets reduced as higher power 
(in 1/t) terms are included. We expect similar effect for gp by allowing for nonzero d 2 and 
f 2 . We determine these coefficients (the 1/t 3 terms) by fitting to the central value of the 
measured ApiT + rate. We obtain —{d 2 + 3 f 2 )/V6 = 1554.6 GeV 6 , which is displayed as the 
second case for Ba in Table ED By assuming d 2 ~ / 2 , we have d 2 —952 GeV 6 , which has 
sign opposite to d \, and is about twice the size of d 2 = f 2 = —478 GeV 6 , the 1/t 3 coefficients 
for the axial vector form factor. 

We show in Fig. d)a) the Apn + decay spectrum. It is interesting that the predicted 
spectra in both d 2 = f 2 = 0 and d 2 — f 2 — —952 GeV 6 cases are close to data. The data 
suggests a curve between these two, which conforms with our expectation that the third, 
1/t 4 term would have same sign as 1/t 2 term. In future as the measured spectrum improved, 
one may in turn use it to extract baryon time-like from factors. 
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FIG. 2: Solid, dashed and dotted lines are for dB(E° £+ vr + )/dm =0 ^qr, dB{E~ S° vr + )/dm„_-^o 
and dB{E~ A7r + )/dm~_-^, respectively, for using (a) the asymptotic gp (cfo = fa = 0), and (b) the 
fitted gp {3,2 = fa = — 952 GeV 6 ). 


While B{Apn + ) is enhanced from the previous results j9LllC| by using our new approach to 
pseudoscalar gp form factor, the enhancement in B(jfa p 7r + ) turns out to be rather mild. This 
can be understood from the relative weight of A vs. S° in Eq. (1A5I) of Appendix A. We expect 
£>(E° pn + ) = 1.6 x 10” 6 , which is within the present Belle limit of B{T,° pn + ) < 3.8 x 10 -6 
at 90% confidence level [u|. Furthermore, it is easy to verify the SU(3) predictions of 
B{T,~ fnr + ) ~ 2 B{Y,° p7r + ) and £>(S°E + 7r + ) ~ 2 B{E~ E°7t + ) given in Table ITT! 

In Fig. [U b) we plot the S°p7r + and T,~ nn + decay spectra. The T,°pn + spectrum is 
close to our previous calculation in Ref. [3|. Since the corresponding SU(3) decomposition 
for these two mode is Dp — Fp, the rates are not sensitive to g ?2 and being zero or finite, 
so long that they are not too different from each other. We show in Fig. [21 the S°S + 7T + , 
5“ S° 7r + and A 7r + decay spectra with c ?2 and fa zero or finite. 

We expect Figs. |T| and [21 to give also the spectra of modes with 7r + replaced by 7r°, but 
with a factor of two reduction in rate from isospin factor. 

In these three-body modes quite often we have a A hyperon produced, which is well 
known to self-analyse its spin upon decay and provides useful information for possible CP- 
and T-violation and chirality studies in B decays Q • Following Ref. [ 27 ]], the angular 
distribution of the cascade B ^ Apn —> ppn decay can be written as 

d 2 T 1 dT 


■[1 + a A {E A ) cos 0], 


(18) 


dE A d cos 6 2 dE A 

where E A is the A energy measured in the B rest frame and 6 is the supplementary angle 
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FIG. 3: (a) a A (E A ), (b) dJ3(Apir + )/dE\ spectrum, where solid (dashed) line is for using the fitted 
(asymptotic) gp of (J 2 = I 2 = — 952 GeV 6 (0). 


between the emitted proton momentum and the B momentum in the A rest frame. We have 
a A {E A ) = V A (E A ) where the A polarization V A (E A ) is given in Appendix B and the 
a A = 0.642 ± 0.013 j2l| is the well-measured A decay asymmetry parameter. 

We show in Fig. 01 the asymmetry a A (E\) and dB(Apir + )/dE A spectrum. The a A (E A ) 
plot is similar to the plot shown in Ref. |27] obtained by using some general arguments. 
The negative a A (E A ) corresponds to a left-handed hclicity dominated A in B decay. It is 
interesting to note that although the decay rate is dominated by the pseudoscalar term, 
we still have polarized A. This can be understood by noting that the ratio of scalar and 
pseudoscalar contributions is roughly given by the averaged fs/dpi which is about 0.1, 
while the polarization V A is roughly given by the averaged —Zfsgp/ifs + 9p) rs./ -2 fs/gp 
which can be as large as —0.6. The sharp turn of a A (E A ) towards much negative value 
for E a > 2.5 GeV is due to the fact that as E A increases, phase space quickly reduces to 
a high m A p region, resulting in the approach of gp to fs and consequently the increase in 
left-handed A polarization. It is well known that the A spin is mainly carried by the s quark 
(as shown in Eq. (A2)) and it is left-handed in the B —> Apn decay (as shown in Eq. 
Therefore, a dominantly left-handed A reflects the V — A nature of the weak interaction 
By comparing Fig. Ha) and Fig. lib), we find — a A ~ 0.2-0.3 for the main portion of Apn 
events. One should be able to check the sign of this asymmetry experimentally in the near 
future. 
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IV. DISCUSSION AND CONCLUSION 


Let us check the 0 3 dependence of the modes considered here. For all modes, By increases 
as we change </> 3 from 60° to 90°; on the other hand, Ba increases for Apn, E N tt and 5“ A n, 


but decreases for 5° E^ tt and S E° tt. However, the variations are at 10 7 order and far 


less significant compared to the Ktt case 


2.8]. Since a§V t bVp s terms dominate, we do not 


expect strong dependence on 0 3 or N c . Similarly, single term dominance implies direct CP 
violation cannot be large, which are found to be within +5% for all modes. 

It is interesting to discuss the implication on ppK and pp tt modes calculated in Ref. jij . 
M . First of all, the changes are in the current-produced parts, whereas these modes contain 
transition parts as well. In particular, the ppTT~ mode is dominated by transition. From 


Eq. 0, we see that /?n(£) is close to its chiral limit form because 
rather weak, and /r^(f) f° r the present work is similar to previous 9, 


lie dependence on gp is 


10|. Therefore, the axial 


vector contributions to ppK and pp tt modes are not affected. The effect of gp only enter 
through the pseudoscalar term. Since the pseudoscalar matrix element for B —> ppK decay, 


{pp |(ss)p|0) 


0 


is Okubo-Zweig-Iizuka (OZI) suppressed, we do not expect much change 


in these modes. On the other hand, for (pp |(dd)p|0) of the ppn mode, by using SU(3) and 


OZI argument as in Ref. 


• 0 . 


corresponds to F P — Dp and is non-negligible. However, this 


0 . 


mode is tree and transition dominant, hence we still do not expect much change in rate 
Note that the transition form factor has a 1/t 3 behavior. For large enough t, the transition 
part is power suppressed. We thus expect to see some 1/t 2 contribution from the new gp 
term, resulting in a slightly broader spectrum than previous jnj . 

In conclusion, we study decay rates and spectra of B —> Apir, E°p7T, E“h7r, 5°E + 7r, 
5“ E° tt and 5“ A tt modes, and the A polarization in this work. By suitably incorporating 
the asymptotic behavior of the baryonic pseudoscalar matrix element, we are able to ob¬ 
tain the ApTT + rate (in part by a fit) and spectrum close to experimental measurements. 


The discrepancy between experimental |Uj and previous theoretical 0 IlD] results is per¬ 
haps resolved. While the ApTT + rate is enhanced from the previous calculation, we expect 
B(T,° pTT + ) = 1.6 x 1CU 6 , which is within the present experimental limit and can be checked 
soon. Although the ApiT + rate is dominated by the pseudoscalar term, we still have A 
polarized giving oa ~ —(0.2-0.3). The impact on ppK due to the present treatment of 
the pseudoscalar form factor is negligible, while we expect a slight broadening of the ppTT~ 
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spectrum. Most of the subtleties in these modes come from the axial and especially the pseu¬ 
doscalar form factors. Information on these form factors may be obtained from studying 
these modes. However, the underlying factorization assumption needs to be checked sepa¬ 
rately. It is interesting that factorization seems to work in ~B° —> D + K~K° and D* + K~K° 
modes, where axial parts are absent and vector parts are known 0 . For these current- 
produced three-body baryonic modes, we expect £>( BB 7r + ) ~ 2H(BB 7r°) as a consequence 
of factorization, which does not depend on the complexity of baryonic form factors. 
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APPENDIX A: ASYMPTOTIC RELATIONS FOR FORM FACTORS f s AND g P 


We follow Ref. 


• 0 


to obtain the asymptotic relations for fs and g P . The wave function 


of an octet baryon can be expressed as 

1 


|B ; T) ~ ^(|B ; UT) + |B ; TTI) + |B ; ITT)), 


(Al) 


i.e. composed of 13-, 12- and 23-symmetric terms, respectively. For B = p, n, E°, A, we 
have 


|p; UT) = 


W y W 3) t(a) _ y| M(1)ti(2)u(3) 


TIT), 


KTIT) = (-b; TIT) with u ^ d), 

|s°;TIT) = 


u(l)d(3) + d(l)«(3) + u{2)d{3) + d(2)u(3) 


V3 

u(l)d(2) + d(l)u{2) 


2VS 


|A;TIT) = 


2\f3 

'd(2)u(3) -u(2)d(3) 


*(3) 


TIT), 




nt>, (A2) 


2 v y 2 
for the corresponding |B ; TIT) parts, while the 12- and 23-symmetric parts can be obtained 
by permutation. To be consistent with the SU(3) decompositions of Table I, our A state has 
an overall negative sign with respect to that of Ref. 
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Following Ref. |l9j, we have 


<B(p)|0|BV)> 

F ± (t) 


u(p ) 


1 + ^ F *(t) + l^2i F -(t) 


2 w 2 
J ±} (0 : B' -»■ B) Fu (t), 


u(p'), 


(A3) 


in the large t limit. Quark mass dependent terms behave like m q /y/\t\ and are neglected. For 
simplicity, we illustrate with the space-like case. Coefficients of F\\ for the O = QlQr, QlQr 
cases are given by 


e+t&^B'-B) = (B; Ut \QW0-, t) - 9(1, i)]|B'; UT> 

+<B; rUI<3[s'(3,t)-9(34)]|B';UT>, 

<i (qLq' R : B' —* B) — 0, 

ettim'L : B' B) = ef : B' B), (A4) 


where Q[</(1(3), |) —> g(l(3),j)] change the parallel spin g'(l(3))| |) part of |B / ;'f||) to a 
g(l(3))| |) part, ft is easy to see that flipping the anti-parallel spin | j) part of |B'; |||) to 
| I) will give a decuplet instead of an octet state. Thus, we need to consider the parallel 
spin case only. By using the above equations, it is straightforward to obtain 


ej"( u L d R :n^p) 
ey" {u l Sr : A —> p) 
et(u L s R : £° —► p) 



(AS) 


By using S, P = Ql^'r ± QrQl an d Eq. (1A4|) . we have ejj z (qq' : B' —► B) = et(qLq' R : B' —> 
B) and e^((pf 5 q' : B' —> B) = ±e^(q L q' R : B' —> B). Hence 

fs = 9 P = eMq L q' R :B'^B)F h (A6) 


in the large t limit. In terms of Dr(p) and Fs(p), we have fs(gp) = Dr(p) + Fs(p), — (-Ds(p) + 
3 Fs(p))/V6, ( D S (p ) —Fs(p))/V2 for B'B = 7ip, A p and £°p cases, respectively. Accordingly, 


Ds = Dp = — Fy , F s = F P = —- F\ 

which implies Eq. (ITU . 


(AT) 
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APPENDIX B: DECAY RATE AND POLARIZATION FORMULA 


For a three-body B —> //BB decay, where h is a pseudoscalar meson and B B is a 
baryon anti-baryon pair, in general the amplitude can be written as 

G f f 


(Bl) 


M(B -> hBB ) = -J= ^Au(pb) h v (PB') + Bu(p B )]/>hl5v(PB') 


The decay rate is given by 


dr = 


+ Cu(p B )v{ps')+ / Du{p B )j 5 v(ps>) 


(£ Ai , 2 \M\ 2 ) dm 2 V2 dm 2 23 , 


(B2) 


(2 7 r) 32 mr B 

where we assign the baryon B as particle 1, the anti-baryon B ; as particle 2 and the meson 
h as particle 3, and Xu 2 ) — ±1 the hclicity of the (anti-)baryon B (B ). 

If the baryon B is in a de fini te helicity state, its spin direction will remain the same in 
either the B meson or its own rest frames. For the baryon B with energy E l (measured in 
the B meson rest frame) the density matrix in the spin (or hclicity) space is given by 


p( E i) = j [! + Pb(Ei)pi -<r] , 


(B3) 


where /3 1 is the unit vector pointing opposite to the direction of the B meson momentum in 
the B baryon rest frame and 

f dm^ Ea 12 (—) Al |A1| 2 


Pb(£i) = 




(B4) 


It is straightforward to obtain: 


s a 1j2 (-) Al \M\ 2 = G 2 F A^Re(AB*)m 1 (2s l -p 3 p2-p3-mls 1 -p 2 ) 
+Re (AT>* — BC*) mim 2 si • p 3 
+Re (AT>* + BC*) (si -p 3 pi ■ p 2 ~ si -pzPi • Ps) 
-Re (C V*) mi si • p 2 \, 


(B5) 


s Ai ,2 l ^| 2 = G 2 f 2 1 \A\ 2 (2pi -p 3 p 2 -p 3 m\pi -P 2 - m 1 m 2 ml) 

+2 Re (AC*) (rri\p 2 ■ p 3 - m 2 pi ■ p 3 ) + |C| 2 (pi • P 2 ~ mim 2 ) 
A —> B, C —> V, m 2 —> —m 2 


(B6) 
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where Si is the hclicity vector of the baryon B (spinor) with Ai = +1. It is easy to check that 
by neglecting rri\ we have rn i .s i —> pi and we obtain Vb{Ei) —> — 1 in the fully left-handed 
chiral case ( A ~ —B and C ~ T>) as expected from Eq. (IB1I) . In general, the polarization 
Vb{E\ ) can be easily evaluated in the B meson rest frame by using 


si = 


(Pb ■ Pi Pi- m\p B ), 


(B7) 


m i\/G °b ■ Pi ) 2 — m\m\ 
where p B is the momentum of the B meson, and the standard technique of expressing p B -pi, 
Pi ■ pj in terms of rrifj. Given these formulas, the task is now reduced to extract the A-T> 
terms for an amplitude of interest. 
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